Accepted to ApJ: August 4, 2010 

Preprint typeset using L£T^]X style cmulatcapj v. 11/10/09 



SIX HIGH-PRECISION TRANSITS OF OGLE-TR-113B 1 

E. R. Adams 2 , M. Lopez-Morales 3 , J. L. Elliot 24 , S. Seager 2 - 4 , D. J. Osip 5 
Accepted to ApJ: August 4, 2010 

ABSTRACT 

We present six new transits of the hot Jupiter OGLE-TR-113b observed with MagIC on the Magellan 
Telescopes between January 2007 and May 2009. We update the system parameters and revise the 
planetary radius to R p = 1.084±0.029 Rj, where the error is dominated by stellar radius uncertainties. 
The new transit midtimes reveal no transit timing variations from a constant ephemcris of greater 
than 13±28 seconds over two years, placing an upper limit of 1-2 on the mass of any perturber 
in a 1:2 or 2:1 mean-motion resonance with OGLE-TR-113b. Combining the new transit epochs with 
five epochs published between 2002 and 2006, we find hints that the orbital period of the planet may 
not be constant, with the best fit indicating a decrease of P = —60 ± 15 milliseconds per year. If 
real, this change in period could result from either a long-period (more than 8 years) timing variation 
due to a massive external perturber, or more intriguingly from the orbital decay of the planet. The 
detection of a changing period is still tentative and requires additional observations, but if confirmed 
it would enable direct tests of tidal stability and dynamical models of close-in planets. 
Subject headings: planetary systems — stars: individual (OGLE-TR-113) 



1. INTRODUCTION 

In the decade since detect i ng th e first transit i ng ex - 
oplanet (jCharbonneau et all 120001: IHenrv et all 120001) . 
over ninety transiting systems have been identified. As 
instrumentation, observations and analysis techniques 
have improved, so has the precision with which param- 
eters of transiting systems can be measured, prompt- 
ing theoretical work to extract more information from 
each light curve. One focus is on how light curves 
are modified by gravitational interactions with other 
planets or tidal interactions with the star. For exam- 
ple, variations in the inclination or duration of succes- 
sive transits would indicate a precessi ng planetary or- 
bit, p otentially from a second planet ()Miralda-Escudel 
120021 ). Periodic midtime var iations could indicate ad- 
ditional planets or moons dHolman fc Murravl 12 005: 
Agol et al.ll2005HHevl & , Gladmanll2007t iFord fc Holman l 
20071 iSimon et all l2007t Kipping) 120091: Kipping et al.l 
2009). A long-term decrease in the orbital period of the 
transiting planet could ind icate orbital decay, a conse- 
quence of tidal dissipation (lSasselovll2003t IPatzold et al.1 
l200i ICarone fc Pa tzold 20071: ILevrard et al J12Q09I). 

In this work we present six new transit light curves 
of the hot Jupiter 0GLE-TR-113b, with a focus on the 
timing effects. In Section 2 we describe the observations, 
data analysis and light curve fitting. Section 3 describes 
our timing analysis and the physical implications of our 
findings. We discuss our results in Section 4. 
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2. OBSERVATIONS 



OGLE-TR-113b was reported bv lUdalski et~aTl (|200l 
as a planet candidate transiting a K-dwarf [1=14.4; 
RA(J2000)=10:52:24.40, Dec(J200 0)=-61:26:48 5l. Its 
pla netary nature was con firmed bv iBouchy et al.l (2004) 
and lKonacki et al.l (|2004[ ). OGLE-TR-113 is an excellent 
target for differential photometry, since it is located in 
a dense region of the sky toward the Galactic plane and 
has numerous nearby, bright comparison stars. We ob- 
served six transits of OGLE-TR-113b as part of a larger 
campaig n to detect t ransit timing variations of OGLE 
planets (|Adamsll2010D . All transits were observed in the 
Sloan j'-band with the dual-CCD instrument MagIC on 
the Magellan Telescopes, located at Las Campanas Ob- 
servatory in Chile. 

Both MagIC CCDs have low readout noise (about 6 e- 
per pixel), small fields-of-view and high-resolution pixels 
(SITe: 142" x 142", or 0."069 per pixel; e2v: 38" x 38", 
or 0."037 per pixel), which both minimizes blends and 
produces stellar images that are spread over many pixels 
(typical FWHM=10-20 pixels depending on seeing and 
binning). This last feature reduces differential pixel re- 
sponse effects and increases the total number of photons 
that can be collected per frame, without requiring defo- 
cusing. The SITe gain was 2.0 e-/ADU, while the e2v 
gain was 2.4 e-/ADU in 2008 and 0.5 e-/ADU in 2009 
(due to engineering changes). The main advantage of 
the e2v is its frame-transfer capability: the readout time 
per frame is only 5s in standard readout mode and 0.003s 
in frame transfer mode, surpassing the 23s readout time 
of the SITe chip and other conventional CCDs. 

We used MaglC-SITe for two transits in January 2007 
and February 2008, henceforth denoted by their UT dates 
as 20070130 and 20080225. MagIC-e2v was used for four 
additional transits between April 2008 and May 2009 
(20080424, 20080514, 20090315, and 20090510). The 
nightly sky conditions ranged from photometric to partly 
cloudy. Therefore, exposure times were adjusted from 10- 
120 seconds to maintain at least 10 6 photons per frame 
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for both the target and one or more comparison stars. 
All transits were observed at air mass < 1.7 to minimize 
differential color effects. The positions of the stars re- 
mained within 5-f pixels each night, while seeing values 
oscillated from 0."4 to 0".7. All transits were observed 
with f x 1 binning, except for 20090510, which was binned 
2x2. Observations lasted from 3.5 to 5 hours to include 
the full transit and out-of-transit baseline. Observations 
were cut short for 20070130 due to clouds. We also trun- 
cated the 20080424 light curve due to a strong systematic 
slope, which we were unable to remove (see Section l2~Tj) . 

Given the utmost importance of accurate timing, we 
took special care to correctly record the time in the 
image headers. The MaglC-SITe and 2008 MagIC-e2v 
times came from the network server, which we verified 
to be synchronized with the observatory's GPS clocks at 
the beginning of each night. The 2009 MagIC-e2v times 
came from an embedded PC104 computer, which receives 
unlabeled GPS pulses every second, and is synchronized 
to the observatory's GPS every night. The intrinsic error 
for all header times is thus much less than a second. 

2.1. Photometry 

All data were overscan-corrected and flatfielded using 
standard IRAF routines^. Simple aperture photometry, 
using the IRAF apphot package, provided optimal re- 
sults, since the target and comparison stars appear well 
isolated on all frames. We ran a wide grid of apertures 
and sky annuli to locate the optimal values, which min- 
imized the out-of-transit photometric dispersion of the 
differential light curves. The comparison stars were iter- 
atively selected to be similar in brightness to the target, 
photometrically stable and un-blended. We examined 
10-20 stars for each transit, using 1, 3, or 10 in the final 
light curves, depending on the night. The best aperture 
radii were 4-24 pixels, while the sky annuli had inner 
radii of 30-120 pixels and 10-pixel widths. A few images 
were discarded because of low counts or saturation. For 
20080424, a 4-pixel aperture (significantly smaller than 
the 15-pixel FWHM radius) was adopted to eliminate a 
slope in the light curve before and during transit, but 
failed to correct the slope after transit. 

We examined the out-of-transit baseline of each light 
curve for systematic trends with respect to air mass, see- 
ing, telescope azi muth, (x,y) pixel loc ation, and time (for 
more details see lAdams et al.l [2010). No trends were 
found in three light curves (20080225, 20080514, and 
20090510), and a fourth (20080424) was not detrended, 
but rather truncated after transit since the slope could 
not be modeled against any variable. We removed a slope 
in telescope azimuth (0.022% deg" 1 ) from the 20070130 
light curve; the transit midtime was unchanged but the 
photometric scatter improved and the depth now agrees 
with the other light curves. We also removed a slope in 
seeing from the 20090315 light curve (—0.16% pixel -1 ), 
observed with variable sky transparency conditions. 

The photometric data for each new transit is excerpted 
in Table [T] All light curves are illustrated in Figure Q] 
binned to 2 minutes for visual comparison, though the 

6 IRAF is distributed by the National Optical Astronomy Ob- 
servatory, operated by the Association of Universities for Research 
in Astronomy, Inc., under cooperative agreement with the National 
Science Foundation 



full data were used in all fits. 

Our analysis also includes four literature light curves: 
20050404 and 20050414 in i?- band (iGillon et al 1 
20050411 in V-b&nd (jPietrukowicz et al 
and 20060318 in if-band 



2006) 

120101 iDi'az et al 



([Snellen fc Coving 120071 ) . The first three light curves 
were kindly provided by the lead author of each paper. 
The original K-band data were lost to a hard-drive crash 
(Snellen; private communication), so we reconstructed 
the unbinned light curve from Figure 3 of their paper. 
We were able to extract the points within a few seconds 
and reconverted the phased data to UTC times with 
zero phase at 2006-03-18 04:56:00 UTC (transit epoch 
1038 from the quoted ephemeris). Finally, we did not 
re-fit the OGLE survey photometry, bu t inste ad used 
the transit epoch from iKonacki et al.1 ([20041 ) in our 
timing analysis (see Section [3]). 

2.2. Light curve fitting 

All lig ht curves in F i gure J were jointly fit us- 
ing the iMandel fc Agoll (120021) a l gorith m, with the 
M CMC method of ICarter fc Winnl ([20091) . as described 
in lAdams et al.l ([2010I ). We assu med quadrati c limb - 
darkening, with initial values from IClaretl ([20001 [2004) , 
using T = 4804 K logg = 4.52, \M /H] = 0, and 
Vmicro — 2 km s _1 ([Santos et al.1 [200 6) . We assumed 
zero limb-dark ening for the flat-bottom ed if-band tran- 
sit, followin g [Snellen fc Covinol ([20071 ) but in contra- 
diction with IClaretl (pOO OH. Throughout we adopted 
M* = 0.78 ± 0.02 Mg,, R* = 0.77 ± 0.02 R Q and 
M p = 1.32±0.19 Mj ([Santos et al.ll2006l ). and zero obliq- 
uity, oblateness and orbital eccentricity for the planet. 
The orbital period was fixed to P = 1.43248 days af- 
ter tests showed little effect on the transit fit results; we 
similarly fixed U2 for each filter. 

The best model fit and corresponding parameter dis- 
tributions were derived from three independent Markov 
chains of 10 6 links each, discarding the first 50,000 
links. We fit values for the planet/star radius ra- 



tio, k = 0.1447 



+0.0006 
-0.0005' 



semimajor axis in stellar radii, 



a/i?» = 6.47_q o9, and orbital inclination, i 



89-0±o.7- 



The fitted transit midtimes are shown in Table [2] We 
derive values for the impact parameter, b = O.lltoogi 
total transit duration T14 = 9647^27 s > planetary radius 
R p = 1.084±0.029i?, 7 (Rj = 71,492 km), and semimajor 
axis, a — 0.02315to;ooogy AU. The errors are discussed in 
Section [231 note that errors on R p and a/i?* incorporate 
stellar radius uncertainties. 

Additionally, each transit was independently fit to 
check for parameter variations over time. Over eight 
years, all values for k, i, a/i?*, and T14 are consistent 
within 1-er of each other and of the joint-fit parameters. 

2.3. Systematic errors 

Light curve systematics caused by atmo- 
spheric/instrumental effects, stellar noise, planetary 
satellites or star spots can affect the correct determi- 
nation of transit parameters. Correlated noise often 

7 The limb-darkening coefficients used are: U\ r = 0.52 ± 0.02, 
u 2 ,r = 0.19 (fixed), iti,v = 0.77 ± 0.04, u 2 ,v = 0.10 (fixed), 
= 0.40 ±0.02, u 2ji / = 0.23 (fixed), u 1}K = (fixed), u 2 , K = 
(fixed). 
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introduces residuals 2-3 times larger than Gaussian 
noise. Such effects have been discussed in detail by 



Panel B shows a new ephemeris fit to all eleven epochs: 



iPont et all (|2006T ) and more recently bv lCarter fc Wirml 
()2009f ) in the specific context of transit midtime 
determination. 

We evaluated how systematics increase the errors es- 
timated from MCMC parameter distributions using two 
meth ods: time-averag e d residuals and re sidual permuta- 
tion (IPont et alj|200€t lSouthworthll2008D . as detailed in 
lAdams et al.l ( 20101 ) . These methods give similar results, 
with errors up to twice those of the MCMC fits, depend- 
ing on the light curve. The errors reported in Table 2 and 
Section POl correspond to the 68.3% credible interval lim- 
its from the MCMC distributions, scaled upwards by the 
excess noise factors for each parameter from the residual 
permutation method. Note that these errors may prove 
in fact to be too conservative in the analysis of transit 
timing var iations in Section 31 

Finally, iSnellen fc CovmT(|2007l ) noted that OGLE- 
TR-113 shows low-amplitude, long-period photometric 
variability, so we have examined whether stellar vari- 
ability could account for correlated features in our light 
curves. For example, the bump before mid-transit in the 
20080514 light curve (Figure[TJ could be interpreted as a 
star spot, but we discard that hypothesis after confirm- 
ing that it coincides with a temporary rapid increase in 
seeing. 

3. TIMING 

To ensure uniform timing analysis, all times in Ta- 
ble [2] are reported in Barycentric Julian Days, using 
the UTC-TT and T T-TDB conver sions (BJD T D B ) af - 
ter lEastman et all (|20Toli 5 l . The iGillon et al.l ([20061 ) 
light curves were supplied with UTC times, which we 
converted to BJDtdb along with our transits. We 
simi larly converted th e recon structed UTC times from 
the ISnellen fc Covinol (120071) light curve. We a dded 
64.184s to the data from iPietrukowicz eFaLl (pOlOh and 
iKonacki et all (|2004D . after confirming that the UTC-TT 
conversion had not been applied. 

Four attempts to fit a transit midtime e phemeris are 
illustr ated in Figure 2. Panel A shows the IGillon et al.l 
(2006) ephemeris, derived using three epoc hs: 20020220 
(OGL E survey), 20050404 and 20050414 (jGillon et al.l 
I2006D . After accounting for the 64.184s UTC-TT con- 
version, the ephemeris is 



T C (N) = 2453464.61740(10)[ J BJD T £, B ] 
+1.43247570(130)^, 



(1) 



where Tc is the predicted transit midtime in BJDtdb, 
the first term is the reference midtime, To, and the second 
term is orbital period, P, times N, the number of elapsed 
transits since To. Although the original data is well fit, 
this ephemeris cannot account for the full data set, giving 
a reduced x 2 =26. 



8 UTC, or Coordinated Universal Time, is typically 
recorded in the image headers. TT, or Terrestrial Time, is 
TT=UTC+32.184s+AT, where AT is the number of leap sec- 
onds since 1972. TDB, or Barycentric Dynamical Time, includes 
millisecond-level relativistic effects. 



T C (N) = 2453464.61762(27)[Bjrj TjDB ] 
+1.43247425(34)^, 



(2) 



with a reduced \ 2 = 1*6. (The errors in Equation 2 
have been re-scaled by Vl.6 = 1.3 to account for the 
deviation of the fit from \ 2 — 1-) Although this equation 
produces a fair fit to the data, it fails to reproduce some 
of the most accurate transit midtimes, e.g. the 20090510 
transit, which deviates by more than 2-a from the fit. 

Panel C shows a third ephemeris fit to only the 2007- 
2009 epochs: 



T C (N) = 2453464.61857(15) [BJD T db] 
+1.43247315(13)^. 



(3) 



This fit has a reduced \ 2 = 0-3 and low errors, but cannot 
account for the 2002-2006 epochs. 

The orbital period in Equation 3 differs by —0.22 ± 
0.11s from the period in Equation 1, with most of the 
error associated with the sparser, lower-precision 2002- 
2005 data. Possible explanations for this period discrep- 
ancy are: (1) the timing errors of the new transits have 
been underestimated — however, the very low reduced 
X 2 of Equation 3 suggests the opposite; (2) the literature 
midtimes contain yet-unaccounted-for systematic errors, 
a possibility which should be investigated, since timing 
errors of hundreds of seconds have been reported before 
(|Adams et alj|201ft IJohnson et al.ll2010D : or (3) the ob- 
served orbital period of OGLE-TR-113b has decreased 
with time. The third explanation, which must be tested 
with additional data, c ould be due to a yet-undetected 
binary star companion (Montalto 2010) with modulation 
period longer than the eight years of observations, or to 
a real, linear decrease in the period due to orbital decay 
of the planet. 

If the change in period is real, it is also rapid. As- 
suming for simplicity a linear rate of change, all 11 tran- 
sit epochs are well fit if the orbital period decreases by 
P = — 60 ± 15 ms yr _1 , with 



T C (N) = 2453464.61873(8) [BJD TDB ] 
(-1.43247426(1 1)N+5P* N(N - l)/2, 



(4) 



where SP = -(2.74 ± 0.66) x 10~ 9 days is the amount 
by which the orbital period changes per orbit, assuming 
P(N) = Pq + SP x N; note P = SP x 365.25/P - The 
reduced x 2 = 0.6 indicates this is a better fit to all epochs 
by a factor of 2.7 than the constant-period fit in Equation 
2. Only the 20050411 epoch deviates by more than 1-a 
(1.7), with residuals of -80 ± 47s. 

As another check on the robustness of this detection, 
we also fit for the average number of periods between 
each pair of midtimes. We find a similar rate, P = — 76 ± 
11 ms yr _1 , with reduced x 2 = 0.7; the fit assuming a 
constant average period has reduced x 2 = 2. 

3.1. Theoretical orbital decay 

We now briefly explore the theoretical implications 
of assuming the observed period decay is real and due 
to orbital decay. Theoretical studies of the closest-in 
planets have until recently focused on OGLE-TR-56b, 
for many years the shortest-period exoplanet known. 
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The estimated timescales of orbital decay vary greatly, 
due to different model assumptions and to uncertainties 
in poorly-constrained parameters. In particular, esti- 
mates of the stellar tidal dissipation factor, Q *, span 
five orders of mag nitude (jPatzold et all 120041 ). Mod- 
els that rely on dissipation by tur bulent v i scosit y find 
lifetimes of a few billion years. ISasselovl (12003ft esti- 
mates a lifetime for OGLE-TR-56b of 0.77 Gyr, imply- 
ing the orbital period changes by P — 2 ms yr" 1 (0.1-5 
ms yr -1 u nder different mode l choic es). iPatzold et al.1 
(2004) and iCarone fc Patzoldl (f2007l ) find a longer life- 
time for OGLE-TR- 56b, 3-3.5 Gyr. A recent study by 
iLevrard et al.l (|2009f ) challenges the assumption that sta- 
ble or unstable tidal equilibrium modes are applicable, 
and finds fast orbital decay rates for most transiting 
planets. Their remaining lifetime for OGLE-TR-56b is 
only 7 Myr, assuming Q* = 10 6 (0.7 Myr to 70 Gyr for 

Q» = 10 5 -10 10 ). 

Adopting the ILevrard et al.l (|2009f ) prediction (their 
equation 5), the lifetime for OGLE-TR-113b is 98 Myr, 
assuming Q* = 10 6 (9.8 Myr to 980 Gyr for Q» = 10 5 - 
10 10 ). Assuming a constant decay rate and that the stel- 
lar Roche lobe radius is reached after the period has de- 
creased by 0.98 days (to Pn oc he = 0.45 d), the fastest 
predicted linear orbital decay rate is 9.8 ms yr^ 1 , if 
Q* = 10 5 . Our tentative period change is six times 
faster, and could be reproduced by the above model with 
Q* = 16000^3000; srx times smaller than the theoretical 
lower estimate. 

3.2. Mass limits on short- period perturbers 

To place upper mass limits on additional objects 
near OGLE-TR-113b, we numerically integrated objects 
wi th periods from 0.3 to 6.3 days, using the approach 
of iSteffen fc Agoll ()2005l ) as described in lAdams et all 
(|2010D . We restricted the analysis to the 6 new tran- 
sits, which can acceptably be fit with a constant period, 
an underlying assumption of the method. In Figure 3 
we show results for perturbers on both initially circular 
(e = 0) and slightly eccentric (e = 0.05) orbits, with 
similar constraints. The strongest constraints are placed 
near the 1:2 and 2:1 mean-motion resonances, where ob- 
jects as small as 1-2 would have been detectable. 
Several additional resonances, including the 5:2 and 3:1, 
also exclude objects larger than 2-20 M®. We therefore 
find no evidence for timing variations caused by close 
companion planets. 

4. DISCUSSION 

We have measured 6 new transits of OGLE-TR-113b, 
which provide the highest-quality data for this planet to 
date. The new transits have timing precisions from 9-28s 
and photometric precisions as good as 0.6-0.7 millimag- 
nitudes in 2 minutes (20080514 and 20090510). The error 
on the refined planetary radius, R p = 1.084 ± 0.029 Rj, 
is dominated by stellar radius uncertainties. 

Observations of OGLE-TR-113b span eight years of 
data (2002-2009) after including five previously pub- 
lished epochs. We find no evidence of short-period timing 
variations, and therefore no sign of other planets in close 
proximity to OGLE-TR-113b. However, there are hints 
that the orbital period calculated from the 2007-2009 
epochs does not agree with orbital period calculated from 



earlier epochs. One explanation for this discrepancy is 
that the period is decreasing by P = —60 ± 15 ms yr" 1 . 

Given the small number of data points, this detection 
is still quite tentative, and systematic effects must be 
carefully considered before any claim is secure. However, 
the fact that the transit parameters are all self-consistent, 
as noted in Section |2~2| suggests that systematics do not 
play a large role. In particular, we find no variation 
greater than 1-er for any of the non-timing parameters 
when each light curve is independently fit. In addition, 
the fits in Equations 2 and 4 have reduced \ 2 values 
significantly smaller than 1.0, indicating that we may be 
overestimating, not underestimating, the timing errors 
in Table 2. 

If the detection is real, one possibility is a wide-period, 
possibly stellar-mass companion producing the pertur- 
ba tions over decade s-long time scales, as recently noted 
by iMontaltol ((2010) . Another possibility is that the or- 
bital period of the planet is decaying. Assuming a linear 
decay rate, the observed rate of period change agrees 
with the theo r etical estimates of planetary lifetime of 
ILevrard et al.l (|2009l ) if the host star has a tidal energy 
dissipation factor Q» = 16000t^^. 

This preliminary hint of orbital period decay needs to 
be verified with observations of the OGLE-TR-113 sys- 
tem, which next becomes observable in January 2011. 
If confirmed, this system has the potential of becom- 
ing a key laboratory to improve current dynamical and 
tidal stability models of close-in planetary systems, and 
to provide observational constraints on Q*. 
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Dave Charbonneau, Dan Fabrycky, Jenny Meyer, Paul 
Schcchter, Brian Taylor, Josh Winn, and the Magellan 
staff for contributions to the observations and analysis; 
Michael Gillon, Pawel Pietrukowicz, Ignas Snellen and 
Andrzej Udalski for helpful communications regarding 
their data; and an anonymous referee for improving this 
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Fig. 1. — Ten transits of OGLE-T R-113b. All transit s, labeled by date and filte r, are binned to 2 minu t es and plotted against orbital 
phase. The top four transits are from Gillon ct al. (2006), Pictrukowicz ct al. (2010), and Snellen & Covino (2007); the rest are new. Solid 
lines show the best joint-fit model, and we report the standard deviation of the residuals. 




Fig. 2. — Observed minus calculated midtimes for OGLE-TR-113b. Literature times are plotted as tri angles and new data as circles; the 
jointly-fit transits are filled symbols. Each panel shows a different ephemeris: (A) constant period from Gillon ct al. (2006), (B) constant 
period using all data, (C) constant period using only 2007-2009 data, (D) linearly-changing period using all data. The curving gray line in 
(A)-(C) shows the changing-period ephemeris from (D). In all panels dashed lines show 1-cr and 3-cr errors on the period; those in (D) also 
incorporate the error on P. 




_i ■ ■ I Li I _J l_ 

1.0 1.5 2.0 3.0 5.0 

Period (days) 



Fig. 3. — Upper mass limit on additional planets with initial eccentricities e c = 0.05 (black) and e c = 0.0 (gray). The central point shows 
OGLE-TR-113b. Near the internal 1:2 mean-motion resonance, obj ects as small as 1-2 wou ld have been detectable. The shaded grey 
region shows the instability region for a 1 Mgj companion, following Barnes & Grccnbcrg (2006). 
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TABLE 1 
Transit Flux Values 3, 



Mid-exposure (UTC) 



Mid-exposure (BJDtdb) 



Flux 



Error 



2454130.599893 

2454130.60128 

2454130.602397 

2454130.603408 

2454130.604375 



2454130.601328 
2454130.602715 
2454130.603832 
2454130.604843 
2454130.60581 



0.9989367 

1.002827 

1.003337 

1.000978 

0.998505 



0.001434 
0.001434 
0.001434 
0.001434 
0.001434 



a This tabic is available in its entirety in a machine-readable form in the online journal. A portion is shown here for guidance regarding its form 
and content. 



TABLE 2 
Transit Midtimes and Residuals 



Transit 


Midtime (BJD TDB ) 


Number 


O-C (s) a 


(7 


O-C (s) b 


a 


20020220 


(2452325.79897 ± 0.00082) c 


-795 


-140 ± 71 


-2.0 


30 ±71 


0.4 


20050404 


2453464.61755 ± 0.00022 





6 ± 19 


-0.3 


0± 19 


0.0 


20050411 


2453471.77900 ±0.00054 


5 


-86 ± 47 


-1.8 


-80 ± 47 


-1.7 


20050414 


2453474.64398 ± 0.00089 


7 


-84 ± 77 


-1.1 


-78 ± 77 


-1.0 


20060318 


2453812.70871 ±0.00046 


7 


-14 ±40 


-0.4 


-28 ± 40 


-0.7 


20070130 


2454130.71844 ± 0.00024 


465 


25 ±21 


1.2 


5 ±21 


0.2 


20080225 


2454521.78377 ±0.00024 


738 


12 ±21 


0.6 


1 ± 21 


0.0 


20080424 


2454580.51523 ± 0.00028 


779 


13 ±24 


0.6 


5 ± 24 


0.2 


20080514 


2454600.56983 ± 0.00011 


793 


11 ± 9 


1.2 


3 ± 9 


0.3 


20090315 


2454905.68672 ± 0.00032 


1006 


-1±28 


-0.0 


13 ±28 


0.5 


20090510 


2454961.55291 ±0.00015 


1045 


-26 ± 13 


-2.0 


-8 ± 13 


-0.6 



a Calculated from constant-period ephemcris Equation [2] 

Calculated from linear ly-changing-period ephemeris Equation [4l 
c Not refit; midtime reported by[Konacki ct al. ( 20 041) . 



